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Abstract 

We provide sufficient conditions on the components of a vector field, which ensure the 
existence of Dulac functions depending on special functions for such vector field. We 
also present some applications and examples in order to illustrate our results. 
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1. Introduction 

Many problems of the qualitative theory of differential equations in the plane are 
related to closed orbits, this fact motivates their study. But deciding whether an arbi¬ 
trary differential equation has periodic orbits or not is a difficult question that remains 
open. 

There are some criteria that allow us to rule out the existence of periodic orbits 
in the plane; between them, we will take particular interest in studying the Bendixson 
Dulac criterion. It is well known that Bendixson-Dulac criterion is a very useful tool 
for investigation of limit cycles (see 0. 0, 0, 0). 

For convenience, we recall the Bendixson-Dulac criterion, see (0, pag. 137). 

Theorem 1.1. (Bendixson-Dulac criterion) Let f) (x \, xi), fi(x \, xf) and h(x\, xz) be 
functions C 1 in a simply connected domain Q C R 2 such that does not 

change sign in Q and vanishes at most on a set of measure zero. Then the system 

[ xi = fi{x\,xf), 

\ X 2 =fi(xi,X 2 ), (X|, x 2 ) e Q, ' 
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does not have periodic orbits in Cl. 

Note that the criterion Bendixson-Dulac also discards existence of poly cycles, mak¬ 
ing it useful in establishing global stability. A function h as in the theorem is called 
a Dulac function. It is well known that Dulac function are an important tool in the 
qualitative study of differential equations, but their determination is a difficult prob¬ 
lem. Our main goal is to give conditions that imply the existence of Dulac functions 
which depends on special forms. We give some consequences and examples to illus¬ 
trate applications of these results. 


2. Results 

Consider the vector field F(x\, xf) = {f\(x \, X 2 ), fi(x \, X 2 )), then the system OJ can 
be rewritten in the form 

x = F(x), x-(xi,X 2 ) e Cl. (2) 

As usual the divergence of the vector field F is defined by 

df\ df 

div(F) - divifufi) = -F— + 

OX 1 OX 2 

We consider C 0 (O, R) the set of continuous functions and define the set 
Fa — 1/ £ C°(£i,R) : / doesn’t change sign and vanishes only on a measure zero set). 

Also for the simply connected region Cl, we introduce the sets 

£&(F) = {h e C‘(aR) : k := > 0, k e Fn) 

OX \ OX 2 

and 

D n (F) = [h £ CHQ.R): k := < 0, k e Fn}- 

OX 1 < 7 X 2 

A Dulac function in the system ([[]) of the Bendixson-Dulac theorem is an element 
in the set 


D n (F) := D + n (F) U D n (F). 

Our results are established with the help of the techniques developed in js}] and 
d, let us recall the following result 


Theorem 2.1. ([9f) If there exist c £ Fn such that h is a solution of the system 
dli dh 

fi-x— + fi-x— = h (c(x 1 , x 2 ) - div(F )), 
dx 1 0 x 2 

with h e Fa, then h is a Dulac function for t(TJ on Cl. 


( 3 ) 
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The next theorem is our first result. In this regard Dulac functions depending on 
special functions 

Theorem 2.2. Let LI be a simply connected open set. Suppose a vector field 
F^f l —+f 2 ^eC 1 (Ll,R 2 ). 

OX 1 OX 2 

If there is c e Ta such that any of the following conditions holds, then Dn(F) £ 0: 

a) .- The function y c-dMF) _ depends on z ■- g\(x\)g 2 (x 2 ) and is continuous; 

JlS 28 i +/ 2 & 1&2 

b) .- The function t) := depends on z. k\(x\) + k 2 (x 2 ) (with k'fxf) = gfxf), 

for i = 1,2) and is continuous; 

c) .- The function cr := c s dn ’[ F \- depends on z. z(x\, x 2 ) and is continuous. 

h dfi + ' 2 jfi 

Proof: We consider the case a), the others are analogous. We seek a Dulac function, 
using the associated equation ([3}. 

First assume that h depends only on z. := gi(xi)g 2 (x 2 ). Thus the equation (0 
reduces to 


f\(xi,x 2 )g 2 (x 2 )g\(xi) 


dh 

dz 


dh 

Mxi,x 2 )gi(xi)g 2 (x 2 )— = h(z)(c(x 1 , x 2 ) ■ 


■ div(F)), 


which is rewritten as 


<9log/; c - div(F) 
dz fg 2 g\ + figig' 2 7 Z ' 

From our hypothesis h := exp ( f y(s)dsj is a solution of the previous equation. Now 

it is easy to verify that h = exp ( j y(s)dsj is indeed a Dulac function. The proof is 
complete. □ 

The following result is a direct consequence of theorem 2.2 and mainly gives some 
particular cases 

Corollary 2.1. Under the conditions of theorem 2.2. If there is c 6 r Fn such that any 
of the following conditions holds, then Dq(F) + 0: 

i) .- The function a ,■ := c dl ^ F) depends only on Xj,for some i e {1,2} and is continuous; 

ii) .- The function f depends on z. x\x 2 and is continuous; 

Hi).- The function 6 := C p+p ) depends on z. := Xi + x 2 and is continuous; 

iv) .- The function e := depends on z. ci-Yi + c 2 x 2 and is continuous; 

v) .- The function k := ffi c( - x i~ x f- d,v( ~F)] depends on z := — and is continuous. 

J Mxi,X 2 )-^f2<.Xl,X 2 ) ^ X 2 

We would like to illustrate the previous results with some examples. 
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Example 1. We consider the classic SIS epidemiological model with disease-induced 
death 

i'l = A - px\ - ax 2 , 

X2 = fi(x 1 - X 2 )x 2 - (a + p + 6 )x2, 

with positive parameters, again we consider K. 2 = {(x\ , X2) € R 2 : X\ > 0, X2 > Oj. 
Denote by F — (f\, f 2 ) the vector field associated to the equation, we get 

div(F) - -p + fixx - 2/3x2 - (a + p + 6), 

taking c(x\ , X2) := —(p + fix 2 ) we can write 

c - div(F) —fix 1 + fix 2 + (a + p + 6) 1 

fi O 2 ) \fix 1 - fix 2 ~(a + p + 5)] x 2 ’ 

which is continuous and depends on z '■= X 2 therefore by i) of corollary 2.1 we get 
£> Rl (F) * 0. 

Example 2. Consider the Lotka-Volterra equations 

xi = X 1 O 1 - kix 1 - b n x 2 ) - hixi = fi(x u x 2 ), 

xi = x 2 (r 2 - k 2 x 2 - biixf) - h 2 x 2 = f 2 (xi,x 2 ), 

since these equations model biological systems in which two species interact, we con¬ 
sider R 2 = {(xi,x 2 ) e R 2 : x\ > 0, x 2 > 0}. Moreover take k\k 2 > 0 with k\,k 2 
are not both zero and we subject each population to their own proportional harvesting 
effort, h\ and h 2 respectively (i.e., the selective harvesting). Denote by F = (f\, f 2 ) the 
vector field associated to the equation, we get 

div(F) = r\ - 2 k\X\ - b\ 2 x 2 - h\ + r 2 - 2k 2 x 2 - b 2 ixi - h 2 , 

taking c(xi,x 2 ) := ~{k\X\ + k 2 x 2 ) we can write 

c - div(F) _ —(r t - k\X\ - bi 2 x 2 - h\ + r 2 - k 2 x 2 - £> 21-^1 - h 2 ) _ -1 

x 2 fi + x\f 2 (xi x 2 ) (r 1 - k\x\ - b\ 2 x 2 - hi + r 2 - k 2 x 2 - b 2 \Xi - h 2 ) x\x 2 ’ 

which is continuous and depends on z ■= X\X 2 therefore by ii) of corollary 2.1 we have 
£>r 

Example 3. Consider the following system of equations 

i'i = (aix 1 + a 2 x 2 )(fii + fi 2 fi' 2 + crix\ p+l ) = f\{x\,x 2 ), 
x 2 = (aixi + a 2 x 2 )(fi 3 +fi 4 x! ” + cr 2 x 2 2 q+l ) - f 2 (x u x 2 ). 

Where a\a 2 > 0 and cr\cr 2 > 0 with <xi, <j 2 are not both zero. If n, m, p and q are non¬ 
negative integers. Denote by F — (fi,f 2 ) the vector field associated to the equation, 
we get 

-div(F) = -affix + fi 2 x J ' 2 + crxx\ p+l ) - cr\{ 2 p + l)xj p (axxi + a 2 x 2 ) 

- a 2 (fii + fi^fif + o- 2 x 2 2 q+l ) - <T 2 (2q + \)x 2 x q {a\xi + a 2 x 2 ), 
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taking c(xi,xi) := (crflp + 1 )x 2 [ p + (T2(2q + l)Xj ? )(ariXi + a2X2) we can write 

c - div(F) 
a\f\ + £*2/2 ’ 

-ai(J3i +/3 2 x n 2 + cr l xl P+1 ) - a 2 (/3i + fax’” + cr 2 x^ +l ) 

{a\X\ + (X2X2) (ui(J3\ +P2X2 + cr\x\ p+X ) + 02(^2 + P^x'" + c^x^ 1 )) 

-1 

a 1 X 1 + 02X2 ’ 

which is continuous and depends on z := a 1X1 + 02X2 therefore by iv) of corollary 2.1 
we have Dr2 (F ) ^ 0 . 


In what follows we present some consequences of Theorem 2.2 and Corollary 2.1 
Let us consider the following system: 


| xi = n(xi)si(x 2 ), 4 

\ x 2 = r 2 (x 1 )s 2 (x 2 ). ' 

We get the next result: 

Corollary 2.2. Let LI c Mr be a simply connected open set. If >'2 si, e Fa cmd r\ 4- 0 
(or r'si e Fa and S 2 ± 0) then system (0 has a Dulac function on LI. 


Proof: We apply i) in Corollary 2.1 □ 

Analogously we can establish the following: 


Corollary 2.3. Let LI be a simply connected open set. Assume that r \ (x \) > 0 (< 0) 
and s' 2 (x 2 ) £ Fa, then system 


x\ = n(xi )r 2 (x 2 ), 

x 2 = Sl(Xl) + s 2 (x 2 ), 


(5) 


has a Dulac function on LI. 


The following produce results on Dulac functions for specific systems 

Proposition 2.1. Let gi : IRC —> R + be continuous functions and a, e IRC for i = 1,2, 
then the planar system 

X\ = Xiffli -xigi(x 2 )), 
x 2 = x 2 (a 2 - x 2 g 2 (x\)), 

supports a Dulac function on M 2 : = {(xi,x 2 ) e R 2 : xi > 0, X 2 > 0}. 
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Proof: First assume a 2 > a 1 and taking c(x 1 , x 2 ) := a 2 — a\ — 2_v 1 ^ 1 (.V 2 J < 0 on E 2 
then the condition i) of corollary 2.1 is written as 

_ c- div(F) 2a 1 - 2xigi(x 2 ) 2 
/1 xi(ai -xigi(x 2 )) xi 

which is continuous and depends on z := x\ therefore by i) of corollary 2.1 we get 
D Rl (F)*<b 

If a\ > a 2 we use a similar argument and get a Dulac function. □ 

Example 4. The following system is a basic model of facultative mutualism M 


X\ 

X 2 


r\x 1 

r 2 x 2 


1 - 

1 - 


xi 

k 1 + b l2 x 2 

x 2 

k 2 + ^2lXl 


/i(xi,x 2 ), 
fi(x i,x 2 ), 


note that satisfies the conditions of the previous proposition, therefore supports a Dulac 
function. 

Similarly we can establish the following result: 

Proposition 2.2. Let gj : R + —> R + be continuous functions and a, 6 R + for i — 1,2, 


xi = xi (a 1 - xigi(x 2 )) - k(x 1 , x 2 ), 
x 2 = x 2 (a 2 - x 2 g 2 (x \)), 

with > 0 and a 2 > a\, then the above system supports a Dulac function on := 
{(xi, x 2 ) e R 2 : x\ > 0, x 2 > 0}. 

Example 5. We consider the following model of two species with harvest Mi 


X'l 

X 2 


rixj 

r 2 x 2 


1 - 

1 - 


Xl 

k\ + b\2X 2 
x 2 

k 2 + b 2 \X\ 


- bx 1 , 


by above proposition, the system supports a Dulac function. 

Proposition 2.3. Let gj : R + —> R + be continuous functions and a t e R + for i =1,2, 
then the planar system 

xi = Xi (gi(Xj) - a,x,), 1 < i + j < 2, 

supports a Dulac function on R 2 := {(xi,x 2 ) el 2 : X! > 0, x 2 > 0}. 

Proof: We verified that the conditions of Corollary 2.1, i) are satisfied, in effect 


-div(F) - —g 1 (x 2 ) + 2a\X\ - g 2 (xi) + 2a 2 x 2 , 
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and taking c = ~(ci\X\ + a 2X1) < 0, then 


c - div(F) 1 

Xlfl+Xlf 2 XiX 2 ’ 

therefore the system supports a Dulac function. □ 

Example 6 . The following system is a model for mutualism (0/ 

X\ = XI [(r, + (r u - n)(l - e~ klX2 )J - fli.ri] = f\(x\,xf), 
x 2 = x 2 \{r 2 + (r 2 2 ~ r 2 )(l - e~ klX ' )) - a 2 x 2 \ = fl(x \, x 2 ), 

where r jt ra, kj, a ,• £ R + are constants and r > r,, i = 1,2. Note that satisfies the 
conditions of the previous proposition, therefore supports a Dulac function. 

Example 7. Gopalsamy /jdj/ had proposed the following model to describe the mutual¬ 
ism mechanism: 


Xl 

X2 


r\X 1 

nx 2 


k\ + a\x 2 

— ;- 1 - Xi 

1 + x 2 

k 2 + 02X1 

— - X2 

1 + Xl 


f\{xi,x 2 ), 

f 2 (xi,x 2 ), 


where r;, kj, a 6 R + are constants and a ,• > kj, i — 1,2. Depending on the nature of kj 
(i = 1, 2), previous system can be classified as facultative, obligate or a combination of 
both. 
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